Long range correlations in the non-equilibrium quantum relaxation of a spin chain 
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We consider the non-stationary quantum relaxation of the Ising spin chain in a transverse field 
of strength h. Starting from a homogeneously magnetized initial state the system approaches a 
stationary state by a process possessing quasi long range correlations in time and space, independent 
of the value of h. In particular the system exhibits aging (or lack of time translational invariance 
on intermediate time scales) although no indications of coarsening are present. 



Non-equilibrium dynamical properties of quantum sys- 
tems have been of interest recently, experimentally and 
theoretically. Measurements on magnetic relaxation at 
low-temperatures show deviations from the classical ex- 
ponential decay which was explained by the effect 
of quantum tunneling. On the theoretical side, among 
others, integrable Q and non-integrable models ||^ were 
studied in the presence of energy or magnetic currents, 
as well as the phenomena of quantum aging in systems 
with long-range Q and short-range interactions 

Here we pose a different question: Consider a quantum 
mechanical interacting many body system described by a 
Hamilton operator H without any coupling to an external 
bath, which means that the system is closed. Suppose the 
system is prepared in a specific state \ipo) at time t = 0, 
which is not an eigenstate of the Hamiltonian H. Then 
we are interested in the natural quantum dynamical evo- 
lution of this state which is described by the Schrodinger 
equation and is formally given by 



|V(t))=exp iV-o) 



(1) 



Obviously the energy E — {ipQ\H\tpo) is conserved. In 
particular we want to study the time evolution of the 
expectation value A{t) of an observable A or the two- 
time correlation function CAB{ti,t2) of two observables 
A and B, defined by 



Ait) = (^oliH(i)IV'o) 

CAB{h,t2) = {^o\{AH{h)BH{t2)}s\^o) 



(2) 



where Anit) — exp{+iHt)Aex-p{—iHt) is the operator 
A in the Hcisenberg picture (with h set to unity) and 
{AB}s — 1/2{AB -f BA) the symmetric product of two 
operators. 

One should emphasize that in such a situation one does 
not expect time translational invariance to hold, which 
would manifest itself in, for instance, A{t) = Aq = const, 
and CABiti^h) = CAsiti —h)- There will be a transient 
regime in which these relations are violated and depend- 
ing on the complexity of the system this non- equilibrium 



regime will extend over the whole time axis, in which we 
would denote it as quantum aging, as it can be observed 
for instance for the universe, which is (most probably) a 
closed system. 

To be concrete we consider the prototype of an inter- 
acting quantum systems, the Ising model in a transverse 
field (TIM) in one dimension defined by the Hamiltonian: 



H = -- 



1 



.1=1 



1=1 



(3) 



where af'^ are spin-1/2 operators on site I. We consider 
initial many body states that are eigenstates either of 
all local af or of all local af operators. We will mostly 
be concerned with fully magnetized initial states, either 
in the x or the z direction, which we denote with |a;) 
and \z), respectively, and which obey aflx) = and 
k) = respectively. 

In passing we note that one obtains the zero temper- 
ature equilibrium situation by choosing the ground state 
of the Hamiltonian (||) as the initial state. This ground 
state has a quantum phase transition at h — I from a 
paramagnetic (h > 1) to a ferromagnetic {h < 1) phase, 
the latter being indicated by long range order in the x- 
component, i.e. a non- vanishing expectation value of . 
Moreover, non-zero ternperature equilibrium relaxation 
has been considered in ||]. 

The expectation values and correlation functions we 
are interested in are those that originate from these spin 
operators af and af . In order to compute them, we 
have to express the Hamiltonian (b|) in terms of fermion 
creation (annihilation) operators rj'^ (rjq) 



(4) 



The energy of modes, eq, q — 1,2, . . . , L are given by the 
solution of the following set of equations 



eq^qil) - -h^q{l)^^qil + l) , 
eq%{l)^-^q{l-l)~h^q{l) , 



(5) 



1 



and we use free boundary conditions which implies for the 
components $g(i + 1) = ^9(0) = 0. The spin-operators 
can then be expressed by the fermion operators as 



af = A1B1A2B2 . 



.Ai.iBi^iAi , 



with 



(6) 



(7) 



and the time-evaluation of the spin operators follows 
from that of the fermion operators: 77^ (i) = e**'^'77+ and 
r]q{t) = e-'*'-'T]g. 

To calculate different non-equilibrium correlation func- 
tions we have developed a systematic method Q in which 
the time-dependent contractions are defined by: 

{AiAk)t = ^cos(e,0*9(0*g(^) , 

{AiBk)t = {BkAi)t = iJ2M^qt)'^qil)-^<,ik) , 

{BiBk)t = J2 cos{e,t)^,il)^,ik) . (8) 

play a central role. For general position of the spin, I = 
0{L/2), one finds simple formulas for the expectation 
values and correlation functions involving operators, 
whereas the calculation of those involving af operators 
is a difficult task and the final result is complicated . 
However, both the surface-spin auto-correlations and the 
end-to-end correlations are given in quite simple form, 
both for the equilibrium |^ and for the non-equilibrium 
case. 

First we study the x-end-to-end correlations defined by 



(9) 



which contain information about the existence or absence 
of magnetic order in the a;-direction. The single time t 
at which the correlations between the two spins are mea- 
sured indicates the age of the system after preparation. 
For the fully ordered initial state |-0o) ~ \x) we obtain 



Cl'^it) = {A,A,)t{BLBL)t + \{AiBL)t 



(10) 



The first term in the r.h.s. of Eq.(lO) is the product 
of surface magnetizations at the two ends of the chain. 
Therefore liuiL^t^oo C^'^ (t) — mi^ and the stationary 
state, starting with \x), has long-range order for /i < 1 
as TOi = 1 ~ h^. Thus the surface order-parameter, mi, 
vanishes continuously at the transition point, he — 1, 
with a non-equilibrium exponent, Z?"*^ — 1. 

The time dependence of the connected correlations 
(generally defined via |2| as CABiti,t2) = CAB{ti-,t2) — 
A{ti)B{t2)) Cl''{t) = \{AiBL)t\^ shows the following 
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FIG. 1. Connected end-to-end correlations Cff (and Cf;^ 
for bottom left) with fixed system size L and field h as a, 
function of the time t calculated with ( |lo|) and (|l^. The left 
column shows data for decreasing field strength h — 1.0, 0.8, 
0.6, the brofcen vertical line is at t = 512 = L. The upper 
right figure shows Cff for different system sizes at H = 1.0, 
the middle right plot shows at /i = 3.0 and the lower right 
figure shows the modulus of for h — 1. Here the broken 
vertical line is at t = 256 — L/2. For the interpretation see 
text. 



features which can be read from fig. ||: 1) They are zero 
for times smaller than a time Th{L) which is equal to the 
system size L for h>l and increases monotonically with 
decreasing h for h <1. 2) Ait = Th{L) a jump occurs to 
a value that decreases algebraically with the system size 
L: 



^max 



(11) 



with a = 2/3 for /i = 1 and a = 1/2 for ft, > 1. 3) 
For t > Th{L) the correlations decay slower than expo- 
nentially, roughly with a stretched exponential. 4) For 
t =^ 3t/i(L) again a sudden jump occurs as for t = Th{L) 
followed by a slightly slower oscillatory decay. 5) This 
pattern is repeated for time t — bTh{L), 7T/i(L), . . ., but 
gets progressively smeared out by oscillations. 

These features can be interpreted as follows: the ele- 
mentary (tunnel) processes of the quantum dynamics of 
the Hamiltonian (^ are spin flips induced by the trans- 
verse field operator a^. In this picture two spins can 
only act coherently and thus give a contribution to the 
connected correlation function if the information about 
such a spin flip processes reaches the two spins simulta- 
neously. Feature 1 tells us that signals generated in the 
center of the system travel with a speed of proportional to 



2 



L/Th{L) to the boundary spins and reaches both simulta- 
neously. At this moment C^^{t) jumps to its maximum 
(see 2). After this, this signal is superposed by other 
more incoherent signals (see 3). However, the strongest 
initial signal is reflected at both boundaries and reaches 
the opposite boundary spins simultaneously again at time 
t = 3Th{L) (see 4), and so on. Of course more and more 
incoherent processes occur in the meantime, giving rise 
to feature 5. 

A similar behavior can be observed for the end-to-end 
correlations when starting with the state jz), which is 



J2 {{AiBk)t{BLAk)t - {AiAk)t{BLBk)t) 



(12) 



The only difference to the behavior of C2'^{t) reported 
above is a) its long time limit vanishes for all values of 
h and b) Th{L), i.e. the earliest time at which the two 
boundary spins are correlated, is only half as big as in the 
previous case. Obviously it is easier to generate and to 
propagate spin flip signals when starting with a z-state. 

Next we study the bulk behavior of the expectation 
values and correlations involving erf operators. We start 
its non-equilibrium expectation value 



■tit) = (^okf (t)IVo) 

^J2{{AiBk)t{BiA(^k))t 

1^ 



(13) 

{AiA,^k))t{BiBk)t) , 



with i{k) — k,{k + 1) for ip = z, (x). We note that 
the equilibrium (i.e. ground state) expectation value, e°, 
corresponds to the energy-density in the two-dimensional 
classical Ising model and we use this terminology also in 
this non-equilibrium situation. For long times the non- 
equilibrium energy-density approaches a finite limit, ef , 
which for a bulk spin is a) for the initial state lipo) ~ \x) 
given by = hjl for /i < 1 and = l/(2/i) for h>\. 
and b) for the initial state |V'o) — by = 1/2 for h < 
1 and = 1 — l/2/i^ for h > 1. Therefore the analogue 
to the specific heat ~ de^/dh is discontinuous at the 
transition point. The relaxation of the energy-density to 
its stationary value is algebraic and follows a i""^/^ low 
for any value of the transverse field, h. At the transition 
point, h = I, we have the analytical results in terms 
of the Bessel-function, J^{x): ef{t) = 1/2 ± Ji(4t)/4t, 
where the -I- (-) sign refer to ip = z(x). 

The two-spin non-equilibrium dynamical and spatial 
correlations involve contributions from different pro- 
cesses described by the contractions ^ and the corre- 
sponding formulas are complicated, therefore they will 
be presented elsewhere Here we report on the basic 
features of the asymptotic behavior of correlations. The 
two-time correlation function {ti <t2), 
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FIG. 2. Connected correlation function C'^'^{r,t) at 
h = 1 given by the expression ( ]l7| ) for different times (i = 2", 
n = 1, 2, 3, . . . from left to right) in a log-log plot. The two 
straight lines indicate the initial dependence of C^'^{r,t) 
for fixed t as well as the dependence of the maximum 
value at r = 2t. 



is non-stationary for intermediate times, ti/{t2 — t\) = 
0(1), which can be read off from our analytical result for 
the connected bulk auto-correlations at /i = 1: 



G 



Z,1p 



{tiM) = Jo (2^2 - 2ii) - i[/(i2 + h)±g{t2 



ti)] 
(15) 



where /(a;) — J2{2x)+Jq{2x), g{x) — J2{2x) — Jq{2x) and 
the -I- (-) sign refer to ip = x{z). Thus we conclude that 
for intermediate times there is aging in the z-component 
auto-correlation function, contrary to what is reported 
in Asymptotically we have limtj^oo Gj^''^(iij ^2) — 
ie^)^, and the connected two-time correlations depends 
only on the time difference, e.g. for h — 1 via ( |T5| ) 
limt.^oo Gf''^(tl, i2) = J§i2[t2 - h]) ~ {Jl(2[t2 - ti])}2. 
For bulk spins this stationary correlation function decays 
algebraically as ^ (^2 — ^i)^^ for any value of h. 

Next we consider the spatial equal-time correlations 

C-'^ir^t) = (Vo|{at./2W'^^./2(0}5|^o) , (16) 

where i = L/2 in a finite system. For long times they 
approach the stationary limit, Muit^oo C^'^ir, t) = {e'^y, 
the same as for the auto-correlation function. For the 
connected correlation function C^''^'{r,t) — C^"^{r,t) — 

^i\^)^t+r^^) '^^'^ derive an analytic expression at the 
transition point, ft, = 1 in the limit L —>■ 00 



J2r(4i) 
~ 1 



(17) 



J2^+l(4i)J2r_l(4t) 



(14) 



which is valid both for \^o) = \x) and \^o) = \z). In fig. |^ 
we show the r-dependence of C^'^{r,t) for various times 
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t. We see that for fixed time t the correlations increase 
proportional to for distances r < t to a maximum 
value C'^ax(0 at r = 2i, which decreases with time pro- 
portional to t~^. For distances larger than r = 2t they 
drop rapidly, faster than exponentially, to zero. 

The latter two features correspond perfectly to what 
we observed also for the z-end-to-end correlations, see 
eq(^2|) : spins that are separated by a distance r can only 
be correlated after the first signal from spin flip processes 
in between them reach simultaneously the two spins, i.e. 
for times t larger than r/2 (for h = 1 and |'0o) = l^}- 
The first feature, that correlations for distances smaller 
than 2t are diminished only algebraically rather than via 
a stretched exponential in the case of end-to-end correla- 
tions, is new and characteristic for bulk spins. For r < 2t 
the correlation function C^''^{r,t) obeys the characteris- 
tic scaling form 



(18) 



with g{x) oa for x ^ 1. The scaling parameter r/t 
appearing in the scaling function g{x) is reminiscent of 
the fact that space and time scales are connected lin- 
early at the critical point in the transverse Ising chain 
since the dynamical exponent is z = 1. Away from the 
critical point we have to evaluate our expressions Q for 
(j^-'^ir^t) numerically for finite but large [L — 512) sys- 
tem sizes. Essentially we observe the same scenario as 
at the critical point, the only difference being that the 
general functional dependency is superposed by strong 
osciUations. Moreover, starting with |^o) = l^^) instead 
of \z) changes the correlations only by a constant factor. 

We collect now our results for the maximum value for 
connected spin-spin correlations since they decay alge- 
braically with various new exponents, we define ourselves 
to /i > 1 since here the time t/j of maximum correlation is 
fixed, whereas for h <1 the value of Th depends on h and 
has to be determined numerically that renders the pre- 
cise determination of the decay exponents difficult. We 
define the ratio a = t/L and amax = Th{L)/L and con- 
sider equal time correlations for fixed values of a. In the 
picture of a propagating front, that separates a region in 
the space-time diagram for the chain in which spins are 
uncorrelated from a region in which they are correlated, 
one observes quasi long range correlations on the front, 
the latter being defined by the ratio t/L — amax- For 
distances smaller than the distance of maximum correla- 
tion or times larger than the correlations decay slower 
than exponential in time, e.g. algebraically for bulk spins 
{C^'^{t,r = fixed) ~ t~^)). When we vary both space 
and time with fixed ratio t/L or t/r we get power laws, 
as long as we stay behind the front (i.e. t > Th)- For 
a > Qfinax we observe again power laws, but with differ- 
ent exponents; they are listed in table 1. 

To conclude we studied a novel type of dynamically 
produced long range correlations in a quantum relaxation 



process in a quantum spin chain. Starting with a homo- 
geneous initial state the quantum mechanical time evo- 
lution according to the Schrodinger equation drives the 
system into a stationary state, which has algebraically 
decaying time-dependent autocorrelations but no criti- 
cal fluctuations. However, during the relaxation process 
spin-spin correlation build up upon arrival of a front of 
coherent signals, which afterwards decay algebraically in 
the bulk. On the front and behind it for fixed ratio of 
space and time scales one observes quasi long range or- 
der. This does not depend on any external parameter like 
the transverse field. This type of algebraic correlation 
needs not to be triggered by some tuning parameter and 
is therefore reminiscent of phenomena in self organized 
criticality |jl^. The scenario we have reported here is a 
result of quantum interference and one may expect that a 
similar one holds for other quantum systems, too. At this 
point one should mention the possibility of coarsening in 
quantum systems as for instance reported in |p^ , which 
is different from the scenario we have reported here. 
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